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I . INTRODUCTION 


In  the  reaction  of  0 atoms  with  OH  two  reactions  are  possible.  One 
is  the  reaction  0+0H  -+•  H+02  (channel  2);  the  other  is  the  exchange  reaction 

0+&  tt*0H+0/  (channel  1).  This  report  will  concentrate  on  the  channel  2 
reaction  but  some  data  on  channel  1 will  also  be  presented. 


There  have  been  a number  of  experimental  measurements  of  the 
thermal  rate  constant  for  the  0+-0H  -*•  H+02  reaction,  all  at  relatively  low 

temperatures  (near  300°K).  The  minimum  thermal  rate  constant  (neglecting 
reaction  complexes,  see  Section  III)  calculated  here  from  trajectories 

agrees  near  300°K  with  several  experimental  results*-3  and  the  recommended 

value3.  However,  this  calculated  rate  constant  is  certainly  low  because 
of  the  long-lived  complex  trajectories  which  the  integrator  of  the 

classical  trajectory  program  (CLASTR^)  cannot  follow  (see  Section  III). 

If  the  decision  of  the  integrator  is  taken  as  correct  for  these  complex 
trajectories,  then  the  rate  constant  calculated  will  be  two  to  three  times 
the  minimum  rate.  If  all  the  complex  trajectories  are  assumed  to  react 
to  form  H+02  (the  worst  possible  case),  then  the  rate  calculated  for 

channel  2 will  be  a little  more  than  three  times  the  minimum  rate. 


Kurzius  and  Boudart  have  formulated  an  empirical  model  to  predict 
the  rate  constant  for  the  0+OH+H+02  reaction  over  a wide  temperature  range 

(see  Fig.  7).  Also  the  reverse  rate  constant  and  equilibrium  constant  data 

harve  been  used  to  predict  the  rate  constant  over  a wide  range  of  temperature 
The  minimum  rate  constant  calculated  here  is  generally  higher  than  these 
7 8 

two  determinations  * (by  a factor  of  four).  The  calculated  minimum  rate 
shows  the  same  flat  temperature  dependence  (Figure  7) . 
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All  the  calculations  have  been  performed  on  a modified  LEPS  surface 
(Surface  I/;)  similar  to  the  surface  (Surface  I)  used  in  the  treatment 

9 

of  the  H+Oj  0H+0  reaction  . Modification  involves  only  a slight  change 

in  the  Sato  parameter  and  addition  of  a constant  to  the  potential  as  a 
whole.  The  potential  contours  of  Surfaces  I and  V'  are  very  similar 
over  all  values  of  the  H-0-0  angle  0. 

The  minimum  thermal  rate  constant  for  the  exchange  reaction 
(channel  1)  is  considerably  less  than  for  the  H+02  channel.  It  varies 

from  one  fifth  to  one  third  the  channel  2 rate. 


II.  POTENTIAL  SURFACE 

The  potential  surface  used  in  these  calculations  is  a modified  LEPS 
surface  similar  to  one  employed  for  the  H+0,  -*■  0H+0  reaction.  The 

10  1 

general  form  of  the  LEPS  surface  is 


V 


V + ^2  + ' (al/2  + “2/2  + °3  2 

- a/  a/  - a/  a/  - c*2'  aj  )H  + d/ 


where 


Q.  D. 

0 ' = -2-  = 

yi  1+A.  4(1+A. ) 


and 


a.  D. 

/ _ i _ i 


i 1+/L  4(1+Ai) 


i = 1,2,3 


(3.al)e■2Si(ri•ri»,-(2.64i)e■6iC^l'rl», 


. . -26.(r.-r.  ) . v -0.(r.-r.  ) 

(l.Jij)  e 11  10-(6*24.)e  11  10 


The  coordinate  system  defining  r^,  r2  and  r3  together  with  certain  other 

variables  is  pictured  in  Fig.  1.  The  zero  of  energy  for  this  surface 
(Surface  V ' ) occurs  for  infinite  separation  in  the  channel  O+OH.  The 
Sato  parameter  is  given  by 


9. 

10. 


0. 53(sin  2Th)2  - 0.25(Cos  Th)2 


A.  Gauss,  Jr.,  J.  Chem.  Phys.  6£,  1689  (1978). 
J.T.  Muckerman,  J.  Chem.  Phys.  56,  2997  (1972). 


2 

♦ 0.2  (V2) 
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where 

b ■ 0.6  X 10'®  cm. 

r2 

Th  » Arcsin  ( — sin  0) 
xa 

- and 

xa  = (r22  * r32/4  ' r2  r3  cos  • 

The  equilibrium  position  for  this  surface  is  very  close  to  that  of 
Surface  I at  r,  = 0.97A,  r-  =>  1.2lA,  0 =«  107°  (H-0-0  angle).  The 
* 11  ^ 

experimental  values  are  r2  =»  0.985A,  r3  = 1.36A,  and  0 = 106°.  The 

variation  of  this  potential  surface  with  the  H-0-0  angle  is  shown  in 
Figure  2.  Also  plotted  here  for  comparison  are  the  ab  initio  results 

of  Gole  and  Hayes1^. 

Potential  energy  contour  plots  are  presented  in  Figures  3,4,  and  5 
for  0 = 180°,  108°  and  80°  respectively.  Note  that  these  contours  have 
the  potential  zero  in  the  separated  H+02  channel  (for  a3  ri  an{* 

r2  -*■  »).  The  zero  of  the  potential  Surface  I77  is  in  the  separated  O+OH 

channel  and  thus  differs  from  the  surface  (called  Surface  I7)*^ 

pictured  in  Figures  3,  4 and  5 by  a constant  (D2e  - D^6) . 

Surface  I had  an  infinity  in  the  derivative  in  the  linear  con- 
figuration. The  change  in  the  angular  (Th)  dependence  for  Surfaces  I7 
and  I77  remedies  this  difficulty. 

The  Surface  I7  (pictured  in  Figures  3,  4 and  5)  has  been  used  for 

dynamics  calculations  of  the  H+C>2  reaction*"*.  Trajectory  results  on 

Surface  I7  closely  match  those  on  Surface  I.  The  system  is  rarely  in 
the  linear  or  near  linear  configuration  even  on  Surface  I7  and  thus 
problems  in  this  configuration  on  Surface  I appear  unimportant  to  the 
trajectory  results. 


11.  J.F.  Ogilvie,  Can.  J.  Spect.,  1_9,  171  (1973). 

12.  J.L.  Gole  and  E.F.  Hayes,  J.  Chem.  Phys.  57^,  360  (1972). 

13.  A.  Gauss,  Jr.,  "The  H+02-*OH+0  Reaction,  A Comparison  of  Potential 

Surfaces",  ARBRL-TR-,  February  1979  (to  be  published). 
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Figure  1 


The  intemuclear  coordinates  r^,  r2  and  r^  and  the  H-0-0 

angle  (0)  are  defined  as  shown.  The  angle  "Th"  and  the 
distance  "x  " are  displayed  on  the  diagram. 

a 
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POTENTIAL  ENERGY 


Figure  2.  Comparison  (at  ^ = 0.96A,  r^  = 1.23A)  of  the  angular 


dependence  of  the  LEPS  potential  energy  Surface  I1 1 with  the 
ab  initio  surfaces  calculated  by  Gole  and  Hayes. 


Figure  3.  Contour  diagram  of  the  LEPS  potential  Surface  I/  for  0=180°. 

Contours  labelled  in  eV.  with  0.5  eV.  spacing.  This  surface 

differs  only  by  an  additive  constant  (D^-D^6)  from  Surface  I7 7 

used  in  this  0 + OH  study. 
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Figure  5. 
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Contour  diagram  (in  eV.)  of  Surface  V for  9=80  . 
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III.  TRAJECTORY  CALCULATIONS  FOR  O+OH 


Trajectory  results  are  summarized  in  Table  I.  All  trajectories  were 
run  for  OH  in  the  ground  state  ( v*0 , J»0) . Table  I displays  the  results 
as  a function  of  relative  translational  energy  (ER).  Three  main  types  of 

trajectory  events  are  recognized:  reactive,  non-reactive,  and  complex. 

The  reactive  events  can  be  subdivided  into  reactive  channel  1 (Rl)  and 
reactive  channel  2 (R2)  events.  Rl  is  the  exchange  reaction,  reaction 
R2  gives  the  H+02  products.  Both  reactive  and  nonreactive  trajectories 

are  back  integrable,  using  the  back  integration  option  of  the  CLASTR 
code.  Both  energy  and  angular  momentum  are  conserved.  Complex  tra- 
jectories are  not  back  integrable,  although  energy  and  angular  momentum 
are  conserved.  In  Table  II  the  complex  trajectories  have  been  further 
subdivided  into  reactive  complexes  and  non-reactive  complexes  according 
to  the  decision  of  the  CLASTR  integrator. 

TABLE  I.  A summary  of  the  reactive  and  complex  trajectories  is  presented. 
The  reactions  are  subdivided  into  channel  1 (Rl)  and  channel 
2 (R2)  reactive  events.  The  total  number  of  trajectories  is 
given  for  each  translational  energy  (ER) . The  OH  is  initially 
in  the  v=0,  J=0  state. 


er 

kcal  \ 

Reactions 

Total 

Complexes 

Total 

Trajs 

mole  I 

5U 

R2 

2 

14 

55 

117 

500 

5 

16 

69 

94 

500 

10 

13 

35 

24 

500 

15 

16 

27 

15 

500 

The  minimum  cross  sections  for  the  channel  O+OH  -+  H+02  are  shown 

in  Figure  6.  Minimum  cross  sections  are  calculated  from  the  R2  events 
only,  complexes  are  assumed  to  be  non- reactive.  The  minimum  cross 
sections  have  here  been  plotted  as  a function  of  translational  energy. 

A simple  straight  line  fits  the  data  well  (Sr  * 0.3  ER  + 6.5).  From 

this  straight  line  fit  the  minimum  specific  rate  constant  (kQO)  has 

been  evaluated  at  four  temperatures  300°K,  400°K,  1600°  and  2500°K 
(See  Table  III).  Evaluation  of  the  specific  rate  constant  at  the  lower 
temperatures  (300  and  400°K)  assumes  that  the  cross  section  has  been 
correctly  extrapolated  into  the  region  of  low  translational  energy 


TABLE  II.  The  complex  trajectories  tabulated  in  Table  I are  here  sub- 
divided into  reactive  and  non-reactive  groups  using  the 
decision  of  the  CLASTR  integrator.  The  complex  reactive 
events  are  further  subdivided  into  channel  1 reactive  events 
(CR1)  and  channel  2 reactive  events  (CR2) . 


er 

/ kcal  \ 
\mole  f 

Total 

Complexes 

Reactive  Complexes  Non-Reactive 

CR1  CR2  Complexes 

2 

117 

11 

76  30 

5 

94 

4 

50  40 

10 

24 

2 

14  8 

15 

15 

2 

3 10 

TABLE  III. 

The  minimum  specific 
for  four  temperatures 

rate 

(300 

constant  for  channel  2 is  shown 
, 400,  1600  and  2500°K) . The 

units  of  the  rate  constant 

here  are  cm  /molecule  sec. 

Minimum 

Temperature  (°K) 

Specific 

Rate  Constant 

3 

(cm  /molecule  sec.) 

300 

(5.4  ±0.7)  X 10'11 

400 

(6.1  ±0.8)  X 10'11 

1600 

(9.3  ±1.2)  X 10"11 

2500 

(9.2  ±1.3)  X 10"11 

(<  2 kcal/mole) . If  rotational  energy  has  little  effect  on  the  specific 
rate  (kQo),  then  it  can  reasonably  be  taken  as  equal  to  the  total  thermal 

rate  constant  (k)  particularly  for  the  lower  temperatures  (<  1000°K) . 
Vibrational  effects  may  have  to  be  considered  at  very  high  temperatures 
but  the  upper  (v  > 0)  vibrational  levels  of  OH  are  not  highly  populated’ 
since  they  are  widely  separated  (>  10  kcal/mole) . Neglecting  rotational 


14 


and  vibrational  effects  the  minimum  specific  rate  constant  kQO  is  taken 


to  be  equivalent  to  the  minimum  overall  thermal  rate  constant  (k)  and  is 
plotted  versus  temperature  in  Figure  7.  Observe  that  the  temperature 
dependence  is  virtually  the  same  as  the  empirical  estimates  of  other 
7 8 

workers  * ; however,  the  minimum  theoretical  rate  is  somewhat  high. 


Another  theoretical  rate  for  the  0+0H 


H+02  reaction  can  be  calcu- 


lated assuming  that  the  decision  of  the  CLASTR  integrator  is  correct 

for  the  complex  trajectories  (Table  III).  The  total  number  of  reactions 

are  then  taken  as  (R2  ♦ CR2) . The  cross  section  curve  for  this  case 

(-BE  2 

S * A R,  A * 17A  , 
r e 

B » 1.22  (kcal/mole)  ‘I  to  the  trajectory  data.  The  total  thermal  rate 
for  this  case  is  shown  in  Figure  7.  The  maximum  possible  rate  (assuming 
all  complexes  are  channel  2 events)  would  be  about  one-third  more  than 
the  rate  assuming  (R2  ♦ CR2)  as  the  correct  number  of  reactions. 


The  CLASTR  program  provides  us  with  some  interesting  information  on 

the  vibrational  and  rotational  excitation  of  the  0.  product  molecule.  At 

kcal  * 

2 relative  translational  energy  most  of  the  energy  goes  into  02 

vibrational  energy  (*"  43%)  with  the  rest  equally  divided  between  rotation 
and  translation  (~  28  to  29%  each).  At  IS  kcal/mole  relative  transla- 
tional energy  most  of  the  energy  goes  into  0 2 rotation  (~  47%)  with 

about  equal  amounts  in  vibration  and  translation  (~  26  to  27%  each). 

The  absolute  energy  in  vibration  stays  nearly  constant  (~  8-9  kcal/mole) 
over  all  translational  energies.  This  vibrational  energy  corresponds 


to  02  excited  to  the  first  or  second  vibrational  level. 


The  exchange  reac.ion  is  much  less  probable  than  the  channel  2 
reaction.  From  the  data  of  Tables  I and  II  cross  sections  may  be  derived. 
The  minimum  cross  sections  for  this  reaction  are  nearly  constant  over  all 


translational  energies;  the  average  is  Sr  > (1.29  t .33)A  . The  cross 


section  points  and  average  are  shown  in  Figure  8.  This  cross  section 
yields  the  minimum  specific  rate  constants  shown  in  Table  IV.  Neglect- 
ing rotational  and  vibrational  effects  these  specific  rates  may  be 
taken  as  the  overall  thermal  rates.  Over  the  temperature  range  300°K  - 
2S00*K  the  exchange  rate  varies  from  one-fifth  to  one-third  the  rate 
in  channel  2 (see  Figure  9). 


(CROSS  SECTION), 


CHANNEL  2 CROSS  SECTIONS 

WITH  CR2  COMPLEXES 

MINIMUM  CROSS  SECTION 


Et  (TRANSLATIONAL  ENERGY),  kcol/mok 

Figure  6.  The  cross  sections  (S  ) for  channel  2 as  a function  of 

r 

translational  energy  (ER)  are  shown.  The  minimum  cross 
section  is  fit  with  a straight  line  Sr  - - 0.3  ER  ♦ 6.5. 

/ -BE» 

The  dashed  line  is  a simple  exponential  fit^Sr»  Ae  , 

A=  1 . B*»0. 122  ^ to  the  data  including  both 

channel  2 reactions  (R2)  and  channel  2 reactive  complexes 
(CR2).  Maximum  impact  parameters  vary  from  3.5  to  4.0A. 
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OTHER , 

KURZIUS  & BOUDART  (1968)7 
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Figure  7 The  rate  constants  (channel  2)  as  a function  of  temperature 

5 7 8 

from  trajectory  calculation  as  well  as  from  other  sources  ’ * 
are  shown.  All  the  rate  constant  curves  show  little  change 
with  temperature  over  the  range  300- 2000® K.  Note  the  units 

3 

cm  /mol.  sec.  on  the  rate  constant. 
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CHANNEL  1 CROSS  SECTIONS 


Z 

KZ 

ttjo 

sfu 

§2 

cc 

£5 

UD 

si 

S5 


H 


/ 

/ 

/ 

hr-® 

/ 


j-V  '(NOI1D3S  SSCRD) 


Figure  8.  The  cross  sections  (ST)  for  channel  1 are  shown  as  a function 

of  translational  energy  (ER) . The  minimum  cross  section  is 

fit  by  the  constant  S = (1.29  ± .33)A2.  The  dashed  line  is 

r -BE 

simple  function  (Sr«A(l+e  R),  A-1.4A2,  B-2.77  (kcal/mole)"1) 

fitting  the  cross  section  data  including  both  channel  1 
reactions  (Rl)  and  channel  1 reactive  complexes  (CR1) . 


(TRANSLATIONAL  ENERGY),  Itcal  /moll 


if 

_ • 
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CHANNEL  1 RATE  CONSTANTS 

WITH  CR1  COMPLEXES 

MINIMUM  RATE  CONSTANT 


103  T"1  /•K'1 


Figure  9.  The  rate  constants  (channel  1)  as  a function  of  temperature 
from  trajectory  calculations  are  pictured.  Note  the  units 

cm  /mol.  sec.  on  the  rate  constant. 


TABLE  IV.  The  minimum  specific  rate  constant  for  channel  1 is  shown  at 
four  temperatures  (300,  400,  1600  and  25008K).  The  units 

of  the  rate  constant  here  are  cm'Vmolecule  sec. 


Temperature  (°K) 

Minimum 

Specific 

Rate  Constant 

3 

(cm  /molecule  sec.) 

300 

(1.1  ± 0.3)  X 10*11 

400 

(2.4  ± 0.6)  X 10"11 

1600 

(3.6  ± 1.0)  X 10"U 

2500 

(4.1  ± 1.1)  X 10"11 

Rates  can  also  be  evaluated  assuming  that  the  decision  of  the 
CLASTR  integrator  is  correct  for  the  complex  trajectories.  The  number 
of  reactions  in  this  case  is  (R1  * CR1). (See  Tables  I and  II.)  The 
cross  sections  for  this  case  are  shown  in  Figure  8.  The  rate  constant 
is  plotted  as  a function  of  temperature  in  Figure  9. 

For  the  channel  1 reaction,  at  2 kcal/mole  relative  translational 
energy  (of  reactants),  only  about  10%  of  the  total  energy  goes  into  vi- 
brational energy  of  the  product  OH,  while  nearly  80%  goes  into  product 
translation.  At  15  kcal/mole  relative  translational  energy  about  30% 
goes  into  vibration  of  the  OH  product  and  some  60%  into  product  trans- 
lation. Product  rotational  energy  is  always  rather  small,  typically 
10  to  15  percent.  Absolute  product  translational  energy  is  nearly 
constant  at  (18  ± 2)  kcal/mole  over  all  relative  translational  energies. 


1 
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IV.  DISCUSSIONS  AND  CONCLUSIONS 

Several  assumptions  have  been  made  in  deriving  estimates  of  the 
thermal  rate  constants.  Among  these  are  the  neglect  of  rotational  and 
vibrational  effects.  These  effects  should  be  analyzed  in  the  future. 
In  particular  the  reaction  0H(v=l)+0  should  be  examined  as  there  has 

14  15 

been  recent  experimental  work  * on  this  reaction. 

14"!  J.E.  Spencer  and  G.P.  Glass,  Intemat.  J.  of  Chem.  Kinetics,  I)(, 
111  (1977). 

15.  J.E.  Spencer,  H.  Endo,  and  G.P.  Glass,  Sixteenth  Symposium  on 
Combustion  (Combustion  Institute,  Pittsburg,  1977)  p.  829. 
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A new  improved  LEPS  potential  surface  for  HO-  has  recently  been 

13  " _ 

derived  . It  uses  two  Sato  parameters  and  fits  the  force  constants  of 

the  HO-  complex  well.  This  surface  has  recently  been  used  for  dynamics 

* 13 

calculations  on  the  H+02  -*■  0H+0  reaction.  A few  preliminary  calcula- 
tions using  this  surface  for  the  O+OH  reaction  indicate  that  it  yields 
cross  sections  similar  to  those  for  the  Surface  l' ' used  in  this  paper. 
However,  many  more  calculations  are  in  order. 

The  low  temperature  rate  constants  (at  300  and  400°K)  are  derived 
by  extrapolating  the  cross  section  curve  to  low  relative  translational 
energies.  Trajectory  calculations  should  be  run  at  these  low  trans- 
lational energies  to  check  the  extrapolated  curve. 

Only  five  hundred  trajectories  were  run  for  each  translational 
energy  because  of  the  large  amounts  of  computer  time  required.  Fifty 
trajectories  at  2 kcal/mole  take  some  700  to  800  sec.  (with  back  inte- 
gration) on  the  CDC  7600.  Even  at  15  kcal/mole,  200  to  300  sec  are 
required  for  fifty  trajectories. 

While  the  minimum  rate  constant  is  in  close  agreement  with  some 
experimental  results  at  300°K  it  is  above  all  quoted  experimental  values. 

1-4 

Four  independent  groups  give  the  rate  constant  at  300  K as 

(5.0  ± 2.0)  X 10'11,  (5.0  ± 2.0)  X lO-11,  (4.3  * 1.3)  X 10'11  and 
-11  3 

(3.3  * 0.5)  X 10  cm  /molecule  sec.;  our  minimum  rate  constant  is 
-11  3 

(5.4  ± 0.7)  X 10  cm  /molecule  sec.  As  pointed  out  earlier  this  cal- 
culated rate  is  definitely  low  because  of  the  complexes,  probably  by  a 
factor  of  two  or  slightly  more.  The  temperature  dependence  of  the 
calculated  rate  constant  closely  follows  the  shape  of  the  Kurzius  and 

Boudart ^ curve  over  all  temperatures. 

All  calculations  have  been  done  on  a single  surface  correlating  with 

3 2 2 

the  ground  state  of  the  reactants  (0(  P)  + OH(  n)),  the  products  (H(  S) 

3 - 2 

♦ 0-  ( Z )),  and  the  HO,  complex  ( A' ' ) . The  complex  has  an  excited 

^ 2 ^ 2 
state  (A')  approximately  17  kcal/mole  above  the  ground  state  ( A/;) 

and  this  excited  state  correlates,  as  does  the  ground  state,  with  the 

3 2 

separated  0(  P)  ♦ 0H(  tt)  reactant  state.  In  the  H+0,  channel  the  excited 

1 2 

state  correlates  with  the  excited  0, ( A)  ♦ H(  S)  state  22.6  kcal/mole 

3 - z 

above  the  ground  state  02(  Z ) ♦ H(  S).  No  attempt  has  been  made  in 

these  calculations  to  account  for  possible  transitions  to  the  excited 
state  although  the  high  rate  of  reaction  (above  all  experimental  values) 
in  channel  2 may  result  from  neglect  of  these  transitions. 
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As  more  data  become  available  on  the  HO 2 surface  through  experiment 

or  ab  initio  calculations  improved  LEPS  semiempirical  surfaces  can  be 
constructed.  In  the  meantime  results  such  as  these  reported  here  help 
to  define  the  effects  of  different  surface  shapes  on  dynamic  calcula- 
tions. 
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